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Abstract 

We study Brcnier's variational models for incompressible Euler equations. These models 
give rise to a relaxation of the Arnold distance in the space of measure-preserving maps 
and, more generally, measure-preserving plans. We analyze the properties of the relaxed 
distance, we show a close link between the Lagrangian and the Eulerian model, and we 
derive necessary and sufficient optimality conditions for minimizers. These conditions take 
into account a modified Lagrangian induced by the pressure field. Moreover, adapting some 
ideas of Shnirelman, we show that, even for non-deterministic final conditions, generalized 
flows can be approximated in energy by flows associated to measure-preserving maps. 

1 Introduction 

The velocity of an incompressible fluid moving inside a region D is mathematically described by 
a time-dependent and divergence-free vector field u{t, x) which is parallel to the boundary dD. 
The Euler equations for incompressible fluids describes the evolution of such velocity field u in 
terms of the pressure field p: 

' dtU+{u-V)u = -Vp m[^,T]xD, 
divM = in[0,T]xi:>, (1.1) 

w • n = on [0, T] X dD. 

Let us assume that u is smooth, so that it produces a unique flow g, given by 

g{t,a) = u{t,g{t,a)), 
g{0,a) = a. 

By the incompressibility condition, we get that at each time t the map g{t, ■) : D ^ D is a 
measure-preserving diffeomorphism of D, that is 



git, ■)#liD = fJ'D, 
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(here and in the sequel f^jj, is the push- forward of a measure ji through a map /, and fiD the 
volume measure of the manifold D). Writing Euler equations in terms of g, we get 



Viewing the space SDiff (D) of measure-preserving diffeomorphisms of D as an infinite-dimensional 
manifold with the metric inherited from the embedding in L^, and with tangent space made by 
the divergence- free vector fields, Arnold interpreted the equation above, and therefore (II. ip . as 
a geodesic equation on SDiff(D) [6]. According to this intepretation, one can look for solutions 
of (|1.2|) by minimizing 



among all paths g{t, •) : [0, T] SDiff(D) with g{0, ■) = f and g{T, ■) = h prescribed (typically, 
by right invariance, / is taken as the identity map i), and the pressure field arises as a Lagrange 
multiplier from the incompressibility constraint (the factor T in front of the integral is just to 
make the functional scale invariant in time). We shall denote by 5{f,h) the Arnold distance 
in SDiff(Z)), whose square is defined by the above-mentioned variational problem in the time 
interval [0, 1]. 

Although in the traditional approach to (11.11) the initial velocity is prescribed, while in the 
minimization of (jl.3p is not, this variational problem has an independent interest and leads 
to deep mathematical questions, namely existence of relaxed solutions, gap phenomena and 
necessary and sufficient optimality conditions, that are investigated in this paper. We also 
remark that no existence result of distributional solutions of (II. ip is known when d > 2 (the 
case d = 2 is different, thanks to the vorticity formulation of (jl.ip ). see |21] . |17| for a discussion 
on this topic and other concepts of weak solutions to (jl.ip . 

On the positive side, Ebin and Marsden proved in [20j that, when D is a smooth compact 
manifold with no boundary, the minimization of (|1.3|) leads to a unique solution, corresponding 
also to a solution to Euler equations, if / and h are sufficienly close in a suitable Sobolev norm. 

On the negative side, Shnirelman proved in [^, [53] that when d > 3 the infimum is not 
attained in general, and that when d = 2 there exists h G SDiff(D) which cannot be connected 
to i by a path with finite action. These "negative" results motivate the study of relaxed versions 
of Arnold's problem. 

The first relaxed version of Arnold's minimization problem was introduced by Brenier in jl2j : 
he considered probability measures rj in 0,{D), the space of continuous paths uj : [0,T] D, 
and minimized the energy 




Vp{t,g{t,a)) {t,a) e[0,T]x D 

a£ D, 
)iS{D) t£[0,T]. 



(1.2) 




(1.3) 




with the constraints 



{eo,eT)#r] = {i,h)#fiD, 



(1.4) 
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(here and in the sequel et{uj) := uj{t) are the evaluation maps at time t). According to Brenier, 
we shall call these ri generalized incompressible flows in [0, T] between i and h. Obviously 
any sufficiently regular path g(t, •) : [0, 1] S{D) induces a generalized incompressible flow 
T] = {^g)^ij,D, where : D ^ Q{D) is given by ^g{x) = g{-,x), but the converse is far from 
being true: the main difference between classical and generalized flows consists in the fact that 
fluid paths starting from different points are allowed to cross at a later time, and fluid paths 
starting from the same point are allowed to split at a later time. This approach is by now quite 
common, see for instance [2] (DiPerna-Lions theory), |TT] (branched optimal transportation), 

Brenier's formulation makes sense not only if h £ SDiff(Z)), but also when h G S{D), 
where S{D) is the space of measure-preserving maps h : D ^ D, not necessarily invertible or 
smooth. In the case D = [0, l]'^, existence of admissible paths with finite action connecting i 
to any h G S{D) was proved in [12], together with the existence of paths with minimal action. 
Furthermore, a consistency result was proved: smooth solutions to (II. ip are optimal even in the 
larger class of the generalized incompressible flows, provided the pressure field p satisfies 

sup sup |V^.p(t,x)| < vr^, (1.5) 

te[o,T] xeD 

and are the unique ones if the inequality is strict. When r) = {^g)^^£) we can recover g{t, •) 
from rj using the identity 

(eo, et)#r] = {i, g{t, ■))#hd, t G [0, T]. 

Brenier found in [12] examples of action-minimizing paths ij (for instance in the unit ball of 
M^, between i and —i) where no such representation is possible. The same examples show that 
the upper bound (II. 5p is sharp. Notice however that {eo,et)#r] is a measure-preserving plan, 
i.e. a probability measure m D x D having both marginals equal to Denoting by T{D) 
the space of measure-preserving plans, it is therefore natural to consider t ^ (eo, 64)^77 as a 
"minimizing geodesic" between i and h in the larger space of measure-preserving plans. Then, 
to be consistent, one has to extend Brenier's minimization problem considering paths connecting 
7, r/ G r(Z)). We define this extension, that reveals to be useful also to connect this model to 
the Eulerian-Lagrangian one in [16] , and to obtain necessary and sufficient optimality conditions 
even when only "deterministic" data i and h are considered (because, as we said, the path might 
be non-deterministic in between). In this presentation of our results, however, to simplify the 
matter as much as possible, we shall consider the case of paths rj between i and h G S{D) only. 
In Section [5] we study the relation between the relaxation 6^ of the Arnold distance, defined 

by 

5,(/i) :=inf |liminf(5(i,/i„) : /i„ G SDiff(L»), / \K - h\'^ dfiD ^ o\ , 
[ Jd J 

and the distance 6{i,h) arising from the minimization of the Lagrangian model. It is not hard 
to show that 5{i,h) < (5^<(/i), and a natural question is whether equality holds, or a gap phe- 
nomenon occurs. In the case D = [0, 1]*^ with c? > 2, an important step forward was obtained by 
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Shnirelman in [24j . who proved that equality holds when h G SDifF(I?); Shnirelman's construc- 
tion provides an approximation (with convergence of the action) of generalized flows connecting 
i to /i by smooth flows still connecting i to h. The main result of this section is the proof that 
no gap phenomenon occurs, still in the case D = [0, l]*^ with d > 2, even when non-deterministic 
final data (i.e. measure-preserving plans) are considered. The proof of this fact is based on 
an auxiliary approximation result, Theorem 15.31 valid in any number of dimensions, which we 
believe of independent interest: it allows to approximate, with convergence of the action, any 
generalized flow r] in [0, 1]*^ by VF^'^ flows (in time) induced by measure-preserving maps g{t, •). 
This fact shows that the "negative" result of Shnirelman on the existence in dimension 2 of non- 
attainable diffeomorphisms is due to the regularity assumption on the path, and it is false if one 
allows for paths in the larger space S{D). The proof of Theorem 15.31 uses some key ideas from 
[24j (in particular the combination of law of large numbers and smoothing of discrete families 
of trajectories), and some ideas coming from the theory of optimal transportation. 

Minimizing generalized paths rj are not unique in general, as shown in [12]; however, Brenier 
proved in [14] that the gradient of the pressure field p, identified by the distributional relation 

Vp{t, x) = —dtVt{x) — div {v (g> Vt{x)) , (1.6) 

is indeed unique. Here Vt{x) is the "effective velocity", defined by {et)^{oj{t)rf) = VtfJ-D, and 
V Vt is the quadratic effective velocity, defined by {et)#{uj{t) <^ t!;(t)?7) = Cg) Vtfio- The proof 
of this fact is based on the so-called dual least action principle: if ij is optimal, we have 

j^tH > -s^Hv) + {P, - (1-7) 

for any measure v in Q{D) such that {eo,eT)#i^ = (i,/i)#/^D and \\p'^ — l\\c^ < 1/2. Here p'^ 
is the (absolutely continuous) density produced by the fiow v, defined by p^{t, ■)pd = {etj^u. 
In this way, the incompressibility constraint can be slightly relaxed and one can work with the 
augmented functional (still minimized by rf) 

I/^^T(i^)-(p,p''-l), 

whose first variation leads to (jl.6p . 

In Theorem 16. 2| still using the key Proposition 2.1 from [Tl], we provide a simpler proof and 
a new interpretation of the dual least action principle. 

A few years later, Brenier introduced in [16] a new relaxed version of Arnold's problem of a 
mixed Eulerian-Lagrangian nature: the idea is to add to the Eulerian variable x a Lagrangian 
one a representing, at least when f = i, the initial position of the particle; then, one minimizes 
a functional of the Eulerian variables (density and velocity), depending also on a. Brenier's 
motivation for looking at the new model was that this formalism allows to show much stronger 
regularity results for the pressure field, namely dxiP are locally finite measures in (0,T) x D. In 
Section [3.31 we describe in detail this new model and, in Section we show that the two models 
are basically equivalent. This result will be used by us to transfer the regularity informations 
on the pressure field up to the Lagrangian model, thus obtaining the validity of (jl.7p for a much 
larger class of generalized flows v, that we call flows with bounded compression. The proof of 
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the equivalence follows by a general principle (Theorem 12.41 borrowed from [3]) that allows to 
move from an Eulerian to a Lagrangian description, lifting solutions to the continuity equation 
to measures in the space of continuous maps. 

In the final section of our paper we look for necessary and sufficient optimality conditions 
for the geodesic problem. These conditions require that the pressure field p is a function and 
not only a distribution: this technical result is achieved in [5], where, by carefully analyzing and 
improving Brenier's difference-quotient argument, we show that dx^p € L^^^ ((0, T); A^ioc(-D)) 
(this implies, by Sobolev embedding, p £ Lf^^(^{0,T); L'l^^:'^ ^\d)^). 

In this final section, although we do not see a serious obstruction to the extension of our 
results to a more general framework, we consider the case of the flat torus T'^ only, and we 
shall denote by /if the canonical measure on the flat torus. We observe that in this case p £ 
Lf^^((0,r);L'^/('^-i)(T'^)) and so, taking into account that the pressure field in (|1.7|) is uniquely 
determined up to additive time-dependent constants, we may assume that fjdP{t, •)dfij = for 
almost all t £ {0,T). 

The first elementary remark is that any integrable function q in (0, T) xT'^ with fj^ Q{t, •) dfij = 
for almost all t £ (0, T) provides us with a null-lagrangian for the geodesic problem, as the 
incompressibility constraint gives 

q{t,uj{t))dtdu{oj) = I / q{t,x) dfij{x) dt = 
nCfi) Jo Jo Jt^ 

for any generalized incompressible flow v. Taking also the constraint (eo, eT)#i^ = {i, h)^^ into 
account, we get 

cT 



£/t{v)=T [ (j \\dj{t)\'^ -q{t,u:)dt] du{uj)> f 
Jn(T<i) \Jo ^ J Jt" 



Cg {x,h{x))dfiTix), 

ln{T'i) \Jo ^ / ^T'' 



where Cg{x, y) is the minimal cost associated with the Lagrangian T ^\uj{t)\'^ — q{t, lo) dt. Since 
this lower bound depends only on h, we obtain that any r] satisfying (II. 4p and concentrated on 
Cg-minimal paths, for some q £ L} , is optimal, and 5 {i,h) = f {i, h) dfij. This is basically 
the argument used by Brenier in [T2] to show the minimality of smooth solutions to , under 
assumption (jl.Sp : indeed, this condition guarantees that solutions of d)(t) = —Vp{t,uj) (i.e. 
stationary paths for the Lagrangian, with q = p) are also minimal. 

We are able to show that basically this condition is necessary and sufficient for optimality 
if the pressure field is globally integrable (see Theorem 16. 12p . However, since no global in time 
regularity result for the pressure field is presently known, we have also been looking for necessary 
and sufficient optimality conditions that don't require the global integrability of the pressure 
field. Using the regularity p £ L^^^ ((0, T); U{D)) for some r > 1, guaranteed in the case D = T'^ 
with r = d/[d — 1) by [5], we show in Theorem 16.81 that any optimal 77 is concentrated on locally 
minimizing paths for the Lagrangian 

^pM := f ^Ht)\^ -Pit,u;)dt (1.8) 
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Since we are going to integrate p along curves, this statement is not invariant under modifications 
of p in negligible sets, and the choice of a specific representative p{t,x) := liminfi^^Q p{t, ■) * 
(peix) in the Lebesgue equivalence class is needed. Moreover, the necessity of pointwise uniform 
estimates on requires the integrability of Mp{t,x), the maximal function of p{t, ■) at x (see 

dSHD). 

In addition, we identify a second necessary (and more hidden) optimality condition. In order 
to state it, let us consider an interval [s, t] C (0, T) and the cost function 

c;'*(x,y) :=inf 1^ ^\ij{T)\^ - p{t,uj) dr : a;(s) = x, u;(t) = y, Mp(r, w) E L^(s, t)| . (1.9) 

(the assumption Mp{t,uj) G L^{s,t) is forced by technical reasons). Recall that, according to 
the theory of optimal transportation, a probability measure A in T'^ x T"' is said to be c-optimal 
if 

/ c{x,y)dX'> / c{x,y)dX 

for any probability measure A' having the same marginals /ii, ^2 of A. We shall also denote 
Wc{ni,fJ-2) the minimal value, i.e. fja^^'' cdX, with A c-optimal. Now, let rj be an optimal 
generalized incompressible flow between i and h; according to the disintegration theorem, we 
can represent r] = J r]^dfj,£){a), with r/^ concentrated on curves starting at a (and ending, since 
our final conditions is deterministic, at h{a)), and consider the plans Aa* = (63,64)^77^. We 
show that 

for all [s,t] C (0,r), A^'* is c^'*-optimal for fij-a.e. a G T'^. (1.10) 

Roughly speaking, this condition tells us that one has not only to move mass from x to y 
achieving Cp*, but also to optimize the distribution of mass between time s and time t. In the 
"deterministic" case when either (eo,es)#r7 or {eo,et)#r] are induced by a transport map g, the 
plan Aa* has Sgi^a) either as first or as second marginal, and therefore it is uniquely determined 
by its marginals (it is indeed the product of them) . This is the reason why condition p.lOp does 
not show up in the deterministic case. 

Finally, we show in Theorem 16.121 that the two conditions are also sufficient, even on general 
manifolds D: if, for some r > 1 and q E L^^^ {{0,T); U' (D)), a generalized incompressible flow 
r/ concentrated on locally minimizing curves for the Lagrangian Cg satisfies 

for all [s,t] C (0,T), A^'* is Cg'*-optimal for ^u/^-a.e. a D, 

then rj is optimal in [0,T], and q is the pressure field. 

These results show a somehow unexpected connection between the variational theory of 
incompressible flows and the theory developed by Bernard-Buffoni [S] of measures in the space 
of action-minimizing curves; in this framework one can fit Mather's theory as well as optimal 
transportation problems on manifolds, with a geometric cost. In our case the only difference 
is that the Lagrangian is possibly nonsmoooth (but hopefully not so bad), and not given a 
priori, but generated by the problem itself. Our approach also yields (see Corollarv l6.13p a new 
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variational characterization of the pressure field, as a maximizer of the family of functionals (for 

[s,t]c{0,T)) 

('?a,7a)d/iT(a), Mq G ([s,t] x T-^) , 
where r]^, 7* are the marginals of Aa*. 

Acknowledgement. We warmly thank Yann Brenier for the many discussions we had together, 
and for the constant support we had from him. We also thank an anonymous referee, who read 
very carefully this paper, for his comments. 

2 Notation and preliminary results 

Measure-theoretic notation. We start by recalling some basic facts in Measure Theory. Let 
X, Y be Polish spaces, i.e. topological spaces whose topology is induced by a complete and 
separable distance. We endow a Polish space X with the corresponding Borel cr-algebra and 
denote by ^{X) (resp. ^^{X), ^{X)) the family of Borel probability (resp. nonnegative 
and finite, real and with finite total variation) measures in X. For A G X and £ ^{X) the 
restriction fi\_A of ^ to A is defined by fi\-A{B) := ^{A n B). We will denote hy i : X ^ X the 
identity map. 

Definition 2.1 (Push-forward) Let /x G ^{X) and let f : X ^ Y be a Borel map. The push- 
forward f^fj, is the measure in Y defined by f^fi{B) = fj,{f~^{B)) for any Borel set B d Y . 
The definition obviously extends, componentwise, to vector-valued measures. 

It is easy to check that has finite total variation as well, and that < An 

elementary approximation by simple functions shows the change of variable formula 

gofdii (2.1) 

for any bounded Borel function (or even either nonnegative or nonpositive, and M-valued, in the 
case 11 G J^+{X)) g:Y 

Definition 2.2 (Narrow convergence and compactness) Narrow (sequential) convergence 
in ^(X) is the convergence induced by the duality with Cb{X), the space of continuous and 
bounded functions in X . By Prokhorov theorem, a family ^ in ^{X) is sequentially relatively 
compact with respect to the narrow convergence if and only if it is tight, i.e. for any e > there 
exists a compact set K C X such that ^{X \ K) < e for any fi G 

In this paper we use only the "easy" implication in Prokhorov theorem, namely that any tight 
family is sequentially relatively compact. It is immediate to check that a sufficient condition 
for tightness of a family of probability measures is the existence of a coercive functional 
^' : A — > [0, +00] (i.e. a functional such that its sublevel sets < i}, i G M"*", are relatively 
compact in A) such that 

/ ^{x)dfj.{x) < 1 \ffi£ ^. 
Jx 
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Lemma 2.3 ([4], Lemma 2.4) Let n G ^(X) and u £ L'^{X;W^). Then, for any Borel map 
f : X , f^{ufi) <^ f^fi and its density v with respect to f^fi satisfies 

/ \v\^df#n < / \u\'^dn. 
Jy Jx 

Furthermore, equality holds if and only if u = v o f fi-a.e. in X . 

Given fi G ^^{X x y), we shall denote by ® A its disintegration via the projection map 
TT{x,y) = x: here A = vr^/U G and x ^ fi^ ^ ^{Y) is a Borel map (i.e. x i-^ f^xi^) is 

Borel for all Borel sets A CY) characterized, up to A-negligible sets, by 

/ f{x,y)dfi{x,y)= [ ( [ f{x,y)dfix{y)] d\{x) (2.2) 

JXxY Jx \Jy J 

for all nonnegative Borel map /. Conversely, any A and any Borel map x^ [i^ G ^(Y) induce 
a probability measure in X x y via (j2.2p . 

Function spaces. We shall denote by ^{D) the space C{[0,T]; D), and by 6<j : [0, T] — > D its 
typical element. The evaluation maps at time t, w i— > to{t), will be denoted by e^. 

If D is a smooth, compact Riemannian manifold without boundary (typically the d-dimensional 
flat torus T'^), we shall denote fin its volume measure, and by do its Riemannian distance, nor- 
malizing the Riemannian metric so that is a probability measure. Although it does not fit 
exactly in this framework, we occasionally consider also the case D = [0,lY, because many 
results have already been obtained in this particular case. 

We shall often consider measures ij £ ^+(0(Z))) such that {et)#r] <^ fio', in this case we 
shall denote hy p'^ : [0, T] x D — > [0, +00] the density, characterized by 

p'^{t,-)fiD:={et)#V, tG[0,T]. 

We denote by SDiff(L') the measure-preserving diffeomorphisms of D, and by S{D) the 
measure-preserving maps in D: 

S{D) ■.= {g:D^D: g#fiD = Pd} ■ (2.3) 

We also set 

S\D) := {g £ S{D) : g is ^^-essentially injective} . (2.4) 

For any g G S^{D) the inverse g~^ is well defined up to -negligible sets, //i^-measurable, and 
g~^ o g = i = g o g~^ p£)-a.e. in D. In particular, if 5 E S'^{D), g~^ G S^{D). 

We shall also denote by T{D) the family of measure-preserving plans, i.e. the probability 
measures in D x D whose first and second marginal are po- 

r{D) := {7 e I^{D X D) : (^i)#7 = /z^, {ti2)#1 = I^d} (2.5) 

(here tti, 112 are the canonical coordinate projections). 
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Recall that SDifF(Z?) C S^{D) C S{D) and that any element g G S{D) canonically induces 
a measure preserving plan 7^, defined by 

7g := {i X g)#HD- 

Furthermore, this correspondence is continuous, as long as convergence in L'^{fi) of the maps g 
and narrow convergence of the plans are considered (see for instance Lemma 2.3 in [3]). Moreover 



{7, : gG S^{D)} = TiD), (2.6) 



SDiff (D)^ ^^''^ = S{D) if D = [0, 1]"^, with d>2 (2.7) 

(the first result is standard, see for example the explicit construction in [181 Theorem 1.4 (i)] in 
the case D = [0,lY, while the second one is proved in |18^ Corollary 1.5]) 

The continuity equation. In the sequel we shall often consider weak solutions fit G ^{D) of 
the continuity equation 

dtfxt + divivtfit) = 0, (2.8) 

where t 1— > /uj is narrowly continuous (this is not restrictive, see for instance Lemma 8.1.2 of 
[3]) and Vt{x) is a suitable velocity field with ||l2(^^) G L^(0, T) (formally, Vt is a section of 
the tangent bundle and \vt\ is computed according to the Riemannian metric). The equation is 
understood in a weak (distributional) sense, by requiring that 

^ / ^it,x)dfit{x)= I dt(l) + {V(l),Vt)dfit inl?'(0,r) 

for any (j) £ ((0)^) ^ D) with bounded first derivatives and support contained in J x D, 
with J (s (0, T). In the case when D C M*^ is compact, we shall consider functions (p € 

((0,T) X R'^), again with support contained in J x R'^, with J d (0,T). 

The following general principle allows to lift solutions of the continuity equation to measures 
in the space of continuous paths. 

Theorem 2.4 (Superposition principle) Assume that either D is a compact subset ofW^, 
or D is a smooth compact Riemannian manifold without boundary, and let fit ■ [0,7"] ^(D) 
be a narrowly continuous solution of the continuity equation ()2.8p for a suitable velocity field 
v{t,x) = Vt{x) satisfying H'l'i |||2(^j) G L^{0,T). Then there exists rj G ^{Q{D)) such that 

(i) fit = {et)#T] for all t G [0,r]; 

(ii) the following energy inequality holds: 

[ [ \dj{t)\'^ dtdriiuj) < [ [ \vt\^dfitdt. 
Jn{D) Jo Jo Jd 

Proof. In the case when D = (and therefore also when C M'^ is closed) this result is 
proved in Theorem 8.2.1 of [3] (see also [7], [25], [10] for related results). In the case when D 
is a smooth, compact Riemannian manifold we recover the same result thanks to an isometric 
embedding in M™, for m large enough. □ 
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3 Variational models for generalized geodesies 



3.1 Arnold's least action problem 

Let f,h^ SDiff (Z?) be given. Following Arnold [6j, we define h) by minimizing the action 

£^t{9) ■.= T [ [ l\g{t,x)\^dfiD{x)dt, 
Jo Jd ^ 

among all smooth curves 

[0,T] 3t^g{t,-) G SDiff(Z)) 

connecting / to h. By time rescaling, b is independent of T . Since right composition with a 
given element g G SDiff(Z)) does not change the action (as it amounts just to a relabelling of 
the initial position with g\ the distance b is right invariant, so it will be often useful to assume, 
in the minimization problem, that / is the identity map. 

The action can also be computed in terms of the velocity field it, defined by x) = 
9{i,y)\y=g-i{t,x)i as 

£/t{u) = tI I ^\u{t,x)\^ dfiD{x)dt. 
Jo Jd 2 

As we mentioned in the introduction, connections between this minimization problem and 
(jl.ip were achieved first by Ebin and Marsden, and then by Brenier: in [12] . |16| he proved 
that if {u,p) is a smooth solution of the Euler equation in [0,T] x D, with D = [0,1]'^, and 
the inequality in (jl.5|) is strict, then the flow g{t, x) of u is the unique solution of Arnold's 
minimization problem with f = i, h = g{T, •). 

By integrating the inequality d'jj{h{x), f{x)) < \g{t,x)\'^ dt one immediately obtains that 
11^ ~ /IIl2(d) — 'v/2(5(/, h); Shnirelman proved in [53] that in the case D = [0, l]*^ with d> 3 the 
Arnold distance is topologically equivalent to the distance: namely, there exist C > 0, a > 
such that 

S{f,g) < C\\f - V/, g G SDiff(D). (3.1) 

Shnirelman also proved in |23j that when d> 3 the infimum is not attained in general and that, 
when d = 2, 6{i, h) need not be finite (i.e., there exist h G SDiff (D) which cannot be connected 
to i by a path with finite action). 

3.2 Brenier's Lagrangian model and its extensions 

In [12], Brenier proposed a relaxed version of the Arnold geodesic problem, and here we present 
more general versions of Brenier's relaxed problem, allowing first for final data in r(D), and 
then for initial and final data in V{D). 

Let 7 G r(L') be given; the class of admissible paths, called by Brenier generalized incom- 
pressible flows, is made by the probability measures rj on Q{D) such that 

(et)#77 = MD ViG[0,r]. 
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Then the action of an admissible rj is defined as 



Jn{D) 

where 

\TL h\uj(t)\'^ dt if is absolutely continuous in [0, Tl 
I +00 otherwise, 

—2 

and 6 (71,7) is defined by minimizing ^xiv) among all generalized incompressible flows r] 
connecting ji to 7, i.e. those satisfying 

(eo,er)#?7 = 7- (3-3) 

Notice that it is not clear, in this purely Lagrangian formulation, how the relaxed distance 
6{7j, 7) between two measure preserving plans might be defined, not even when r] and 7 are 
induced by maps g, h. Only when g £ S^{D) we might use the right invariance and define 

%3,7/i) :=%i,7ho<;-i)- 

These remarks led us to the following more general problem: let us denote 

n{D) := n{D) X D, 

whose typical element will be denoted by {u>,a), and let us denote by ttd : 0,{D) D the 
canonical projection. We consider probability measures r} in ^[D) having as second marginal, 
i.e. (7r£))^77 = they can be canonically represented as r]^ ® fj^D, where r]^ £ ^{Q{D)). The 
incompressibility constraint now becomes 

/ {et)#VadMa) =t^D Vt E [0,r], (3.4) 

JD 

or equivalently {et)^T) = for all t, if we consider ej as a map defined on n{D). Given initial 
and final data ij = r]a fJ-o, 7 = 7a'S>>/UD G ^(Z)), the constraint (|3.3p now becomes 

ieo,iTD)#r] = Va'S' fJ-D, {eT,T^D)#ri = la® P'D- (3.5) 
Equivalently, in terms of 77^ we can write 

(eo)#?7a = ^a, {eT)#Va = la- (3.6) 

_2 

Then, we define 6 (?y, 7) by minimizing the action 

s^t{'-^) dri{uj, a) 

n(D) 

among all generalized incompressible flows r/ (according to (j3.4p ) connecting to 7 (according 
____ 2 

to (j3.5|) or (|3.6|) 1. Notice that S is independent of T, because the action is scaling invariant; 
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so we can use any interval [a, b] in place of [0, T] to define S, and in this case we shall talk of 
generalized flow between rj and 7 in [a, b] (this extension will play a role in Remark 13.21 below) . 

When rja = da (i.e. r/ = 7i), ()3.6p tells us that almost all trajectories of ij^ start from a: 
then Jj^ T]^ d^£,[a) provides us with a solution of Brenier's original model with the same action, 
connecting 7i to 7. Conversely, any solution v of this model can be written as Uadfiu, with 
Va concentrated on the curves starting at a, and fa® ^i-n provides us with an admissible path 
for our generalized problem, connecting 7^ to 7, with the same action. 

Let us now analyze the properties of {T{D), 5)] the fact that this is a metric space and even 
a length space (i.e. any two points can be joined by a geodesic with length equal to the distance) 
follows by the basic operations reparameterization, restriction and concatenation of generalized 
flows, that we are now going to describe. 

Remark 3.1 (Repameterization) Let x '■ [0)^] ~^ [0)^] t>e a map with x > 0, x(0) = 
and x(^) = ^- Then, right composition of u! with x induces a transformation rj 1— > x#'? 
between generalized incompressible flows that preserves the initial and final conditionl. As a 
consequence, if rj is optimal the functional x ^ -^Tix^v) attains its minimum when xit) = t- 
Changing variables we obtain 

s^T{x#r]) = T f xHt) f \\u\\x{t))dr]{uj,a)dt = T f / {s) dr^{uj , a) ds 

Jo Jn(D) ^ Jo 9[s) Jn{D) ^ 

with g = x~^- Therefore, choosing g{s) = s + £(p{s), with cj) G C^(0,T), the first variation gives 

[ ( [ \dj\'^is)dri{u;,a)) <j){s)ds = 0. 
Jo \Jh{D) J 

This proves that s ^ /^^(d) {^{"^(s) dr]{uj,a) is equivalent to a constant. We shall call the square 
root of this quantity speed of rj. 

Remark 3.2 (Restriction and concatenation) Let [s,t] C [0,T] and let r<i^f : C{[0,T]; D) — > 
C{[s,t];D) be the restriction map. It is immediate to check that, for any generalized incom- 
pressible flow T] = T}^® fiD in [0, T] between r/ and 7, the measure {rs^t)#'n is a generalized 
incompressible flow in [s,t] between r/^ := {es)#r]^ (g) fio and 7^ := {et)#r]^ fi£), with action 
equal to 

(t-s) / -\uj{t)\^ dTdr]{uj,a). 

Jn{D) Js ^ 

Let s < I < t and let r] = fj,g^(dfJ^D, = Ua(^^XD be generalized incompressible flows, respectively 
defined in [s,/] and and joining 77 to 7 and 7 to 6. Then, writing 7^ = {ei)^r]^ = [ei)^Ua, 
we can disintegrate both rj^ and Va with respect to 7^ to obtain 

Va = I Va,. djaix) G ^{C{[s, l];D)), Ua= [ u^,, dja{x) S ^{C{[1, t];D)), 
Jd J d 

with J., Ua^x Concentrated on the curves u with uj{l) = x. We can then consider the image Xx,a, 
via the concatenation of paths (from the product of C{[s,l];D) and C{[l,t\]D) to C([s,t];D)), 
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of the product measure rjax ^ ^a,x to obtain a probability measure in C{[s,t\;D) concentrated 
on paths passing through x at time /. Eventually, setting 




we obtain a generalized incompressible flow in [s, t] joining rj to 9 with action given by 

where -^[3,1] iv) is the action of 77 in [s, /] and -^[i^t] (^) is action of v in [I, t] (strictly speaking, 
the action of their restrictions). 

A simple consequence of the previous remarks is that 5 is a distance in T{D) (it suffices to 
concatenate flows with unit speed); in addition, the restriction of an optimal incompressible flow 
77 = 77^ (g) between rja® fJ-D and 7a (8) //d to an interval [s, t] is still an optimal incompressible 
flow in [s, t] between the plans {es)^r]^ ® fiD and {et)#r]^ (g) This property will be useful in 
Section [6l 

Another important property of 5 that will be useful in Section [6] is its lower semicontinuity 
with respect to the narrow convergence, that we are going to prove in the next theorem. Another 
non-trivial fact is the existence of at least one generalized incompressible flow with finite action. 
In [12^ Section 4] Brenier proved the existence of such a fiow in the case D = T'^. Then m 
[241 Section 2], using a (non-injective) Lipschitz measure-preserving map from T*^ to [0, 1]"^, 
Shnirelman produced a fiow with finite action also in this case (see also Section 3]). In the 
next theorem we will show how to construct a flow with finite action in a compact subset D 
whenever flows with flnite action can be built in D' and a possibly non-injective, Lipschitz and 
measure-preserving map f : D' ^ D exists. 

Theorem 3.3 Assume that D CR"^ is a compact set. Then the infimum in the definition of 
S{rj,^) is achieved, 

(rj,^) 5{r],'j) is narrowly lower semicontinuous (3-7) 

and 

5(7i,7h) < S{i,h) V/i G SDiff(i?). (3.8) 

Furthermore, sup 6{rj, 7) < \/d when either D = [0, 1]*^ or D = T'^ and, more generally, 
r,,7er(D) 

sup 5z)(7i,7) < Lip(/) sup 6D'{ji,j') 
7er(D) 7'Gr(D') 

whenever a Lipschitz measure-preserving map f : D' ^ D exists. 

Proof. The inequality 5(7i,7/i) < 5{i,h) simply follows by the fact that any smooth flow g 
induces a generalized one, with the same action, by the formula rj = where <I> : Z) — > r2(D) 
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is the map x i— > {g{-,x),x). Assuming that some generahzed incompressible flow with a finite 
action between rj and 7 exists, the existence of an optimal one follows by the narrow lower 
semicontinuity of 77 1— > s^Tiff) (because lo i— > s^t{^) is lower semicontinuous in ^{D)) and by 
the tightness of minimizing sequences (because £/t{^) is coercive in Vl{D)^ by the Ascoli-Arzela 
theorem). A similar argument also proves the lower semicontinuity of (??, 7) ^ 5{ri,j), as the 
conditions (j3.4p . (j3.5p are stable under narrow convergence (of rj and rj, 7). 

When either D = [0, l]'^ or D = T"^, it follows by the explicit construction in |12j . [24j 
that 5(7i,7/i) < for all h G S{D); by right invariance (see Proposition 13.41 below) the same 
estimate holds for (5(7j,7ft,) with / G S'^{D); by density and lower semicontinuity it extends to 
^(?7,7), with ?7, 7 e r(D). 

Let / : D' — > D be a Lipschitz measure-preserving map and h € S{D); we claim that it 
suffices to show the existence of 7' G ^{D') such that (/ x /)#7' = (i x h)^^D- Indeed, if this 
is proved, since / naturally induces by left composition a map F from Q,[D') to 0(1?) given by 
(w(t),a) I— > (/(a;(t)), a), then to any G Q{D') connecting i to 7' we can associate F^rj, which 
will be a generalized incompressible flow connecting i to h. By the trivial estimate 

^T{F#'n) < Lip2(/)i2/r(r?), 

one obtains Sci{^i,h) < Lip(/)(5£)/(7^, 7'). By density and lower semicontinuity we get the 
estimate on 60(71,^) for all 7 G r(L'). 

Thus, to conclude the proof, we have to construct 7'. Let us consider the disintegration of 
fiD' induced by the map /, that is 



f^D' = / f^ydfiniy) (3.9) 



D 



where, for /i/j-a.e. y, fiy is a probability measure in D' concentrated on the compact set / ^{y). 
We now define 7' as 

7 := / flyX fih{y)dfJ'D{y)- 



D 



Clearly the first marginal of 7' is fio'] since h G S{D), changing variables in ()3.9p one has 
Hd' = Jd l^h{y) dfJ-oiy), and so also the second marginal of 7' is fiD- Let us now prove that 
(/ X /)#7' = (i X /i)#/iD: for any (/> G Cfe(D x Z?) we have 

cPiv, y') d{f X /)#7'(y, y') = [ <^(/(x), /(x')) dj\x, x') 
DxD Jd'xD' 




4>{f{x), f{x')) dfiy{x) dfi^y) {x') duoiy) 

D JD'xD' 

4>{y^Hy))dfiD{y), 



D 



where in the last equality we used that fiy is concentrated on / ^(y) and fJ-hiy) is concentrated 
on f~^{h{y)) for ^zj-a.e. y. □ 
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By (13. ip . ()3.8p and the narrow lower semicontinuity of 6{i, ■) we get 

< C||/i-i||22(B) if /i G 5(Z)), i:* = [0,1]"^, d> 3. (3.10) 

We conclude this section by pointing out some additional properties of the metric space 
{T{D),6). 

Proposition 3.4 (T{D),6) is a complete metric space, whose convergence implies narrow con- 
vergence. Furthermore, the distance 6 is right invariant under the action of S^{D) on T{D). 
Finally, 6 -convergence is strictly stronger than narrow convergence and, as a consequence, 
(r(D),(5) is not compact. 

Proof. We will prove that 5{rj,"f) > W2(?7, 7), where W2 is the quadratic Wasserstein dis- 
tance in ^3^{D X D) (with the quadratic cost c((xi, X2), (yi, 2/2)) = c^dC^^i, yi)/2 + (i|)(x2, y2)/2); 
as this distance metrizes the narrow convergence, this will give the implication between 5- 
convergence and narrow convergence. In order to show the inequality 6{r],^) > W2{r],^) we 
consider an optimal flow r]^ defined in [0, 1]; then, denoting by uJa G ^{D) the constant 
path identically equal to a, and by Va £ =^(C([0, 1];D x D)) the measure 77^ x 5i^^, the measure 

u := Ua dfj,D{a) G =^(C([0, 1]; -D x D)) provides a "dynamical transference plan" connecting rj 

2 

to 7 (i.e. (eo)#^' = ??, {ei)^u = 7, see (271 Chapter 7]) whose action is 6 (r/, 7); since the action 

of any dynamical transference plan bounds from above (^)7)) the inequality is achieved. 

The completeness of {T{D),5) is a consequence of the inequality 6 > W2 (so that Cauchy 
sequences in this space are Cauchy sequences for the Wasserstein distance), the completeness 
of the Wasserstein spaces of probability measures and the narrow lower semicontinuity of 5: we 
leave the details of the simple proof to the reader. 

The right invariance of 5 simply follows by the fact that rjok = r]ii(^g^-^ ®^d, 70/1 = 7/^(3) (^fJ^D, 
so that 

6{r] o /i, 7 o /i) < 5{ri, 7), 

because we can apply the same transformation to any admissible flow rj^^i^ fi^ connecting r] to 
7, producing an admissible flow r]^^-^ ® l^D between rj o h and 70/1 with the same action. If 
h G S'^{D) the inequality can be reversed, using h~^. 

Now, let us prove the last part of the statement. We first show that 

I [ dl{f,h)d^iD<6\lf^h) yf,heSiD). (3.11) 

Indeed, considering again an optimal flow rj^^® ^lo, for fi£)-a.e. a £ D we have 

ldl{f{a),h{a)) = VF|('5/(a),5,(,)) <T [ f \\u{t)\^ dt drj^iu), 
^ Jn{D) Jo ^ 

and we need only to integrate this inequality with respect to a. From (j3.1ip we obtain that S'(-D) 
is a closed subset of r(Z)), relative to the distance b. In particular, considering for instance a 
sequence (3^1) C S{I)) narrowly converging to 7 G r(Z)) \ S{D), whose existence is ensured by 
(j2.6p . one proves that the two topologies are not equivalent and the space is not compact. □ 
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Combining right invariance with ()3.10p . we obtain 

Xlgnh)=~6{li,lhog-^)<C\\g-h\\l,^j,^ WheSiD), geS\D) (3.12) 

a D = [0,1]^ with d > 3. By the density of S^{D) in S{D) in the norm and the lower 
semicontinuity of 6, this inequahty stih holds when g G S{D). 

3.3 Brenier's Eulerian-Lagrangian model 

In [16], Brenier proposed a second possible relaxation of Arnold's problem, motivated by the 
fact that this second relaxation allows for a much more precise description of the pressure field, 
compared to the Lagrangian model (see Section [6]). 

Still denoting by = ?7a Hd £ r(Z)), 7 = 7a fJ-D G r(Z)) the initial and final plan, 
respectively, the idea is to add to the Eulerian variable x a Lagrangian one a (which, in the 
case rj = 'ji, simply labels the position of the particle at time 0) and to consider the family of 
distributional solutions, indexed by a G D, of the continuity equation 

dtct,a + dw{vt,act,a) = in V'{{0, T)xD), for ^D-a.e. a, (3.13) 

with the initial and final conditions 

co,a = ila, CT,a = 7a, for no-a.e. a. (3.14) 

Notice that minimization of the kinetic energy Jj^ \ vt,a\'^ dct^a dt among all possible solutions 
of the continuity equation would give, according to [8], the optimal transport problem between 
r]a and 7a (for instance, a path of Dirac masses on a geodesic connecting g{a) to h{a) if rja = 6g(^a) ) 
7a = (5/i(a))- Here, instead, by averaging with respect to a we minimize the mean kinetic energy 

I \Vt,a\'^ dct^adtdfiD{a) 
D 

with the only global constraint between the family {ct^a} given by the incompressibility of the 
flow: 

[ ct.adfiD{a)=fiD VtG[0,r]. (3.15) 
Jd 

It is useful to rewrite this minimization problem in terms of the the global measure c in [0, T] x 
D X D and the measures ct in D x D 

C ■■= Ct,a ® (-Sf^ X lln), Ct := Q,a fJ-D 

(from whom ct,a can obviously be recovered by disintegration), and the velocity field v{t, x, a) := 
Vt,a{x)- the action becomes 

£/t{c,v):=T I I -\v{t,x,a)\'^ dc{t,x,a), 
Jo JdxD 2 
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while (I3.13P is easily seen to be equivalent to 



— / 4>{x,a)dct{x,a) = / {Vxcl){x,a),v{t,x,a)) dct{x,a) (3.16) 
JdxD JDxD 

for all (j) £ Cb{D X D) with a bounded gradient with respect to the x variable. 

Thus, we can minimize the action on the class of couples measures- velocity fields (c, v) that 
satisfy ()3.16p and ()3.15p . with the endpoint condition ()3.14p . The existence of a minimum in 
this class can be proved by standard compactness and lower semicontinuity arguments (see [16] 
for details). This minimization problem leads to a squared distance between r] and 7, that we 
shall still denote by 6 (t/, 7). Our notation is justified by the essential equivalence of the two 
models, proved in the next section. 



4 Equivalence of the two relaxed models 

In this section we show that the Lagrangian model is equivalent to the Eulerian-Lagrangian 
one, in the sense that minimal values are the same, and there is a way (not canonical, in one 
direction) to pass from minimizers of one problem to minimizers of the other one. 

Theorem 4.1 With the notations of Sections and \3.3l 

min^/Tiv) — niin=e^(c, i') 

r? {c,v) 

for any rj, j £ ^(D). More precisely, any minimizer T] of the Lagrangian model connecting r] to 
7 induces in a canonical way a minimizer (c, v) of the Eulerian-Lagrangian one, and satisfies 
for -a.e. t G [0,T] the condition 

uj{t) = vt,a{et{(^)) forrj-a.e. {uj,a). (4.1) 

Proof. Up to an isometric embedding, we shall assume that D C M*" isometrically (this is 
needed to apply Lemma [2^ . If rj = rj^ fj, £ ^{^1{D)) is a generalized incompressible flow, 
we denote hy D' C D a Borel set of full measure such that ^/xiVa) < ^ ^'^^ a £ D'. For any 
a £ D' we define 

c^a := {et)#Va, ^7,a = {'^it)Va) ■ 

Notice that m^^ is well defined for ^^-a.e. t, and absolutely continuous with respect to c^^, 
thanks to Lemma [2.3| moreover, denoting by v^^^ the density of with respect to cj^^, by the 
same lemma we have 

/ \v?,a\'dcl,< I \u{t)\''dr^S-^), (4.2) 
Jd Jn{D) 

with equality only if uj{t) = v'^^{et{uj)) for r/^-a.e. uj. Then, we define the global measure and 
velocity by 

c"' := cj^a ® i-^^ X IJ-d), v'^{t,x,a) = v^{x,a) := v'l',^{x). 
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It is easy to check that {c^,v'^) is admissible: indeed, writing rj = rja fi, j = 'ja ^J-d, the 
conditions (eo)#??„ = ??a a nd {eT)#Va = la yield c^^„ = r?a and c^^^ = 7a (fo r ^D-a.e. a). 

This proves that ()3.14p is fulfilled; the incompressibility constraint ()3.15p simply comes from 
p.4p . Finally, we check (j3.13p for a ^ D'\ this is equivalent, recalling the definition of Vt^a, to 



'dt 



<l,{x)dcl,{x)= / (V</>,m,''J, (4.3) 



which in turn corresponds to 



d 
It 



<t>{u{t))drjaiu^)= / {V^{uj{t)),d;{t))drj,{uj). (4.4) 
n(D) Jn{D) 



This last identity is a direct consequence of an exchange of differentiation and integral. 

By integrating (j4.2p in time and with respect to a we obtain that iS^(c^,i)^) < s^Ti'n)^ 
equality holds only if ()4.ip holds. 

So, in order to conclude the proof, it remains to find, given a couple measure- velocity field 
(c, 1)) with finite action that satisfies ()3.13p . (I3.14p and ()3.15p . an admissible generalized incom- 
pressible flow r/ with ^/xiv) — (£,!>). By applying Theorem 12.41 to the family of solutions of 
the continuity equations (|3.13p . we obtain probability measures with {et)^r]^ = ct^a and 

/ \uj{t)\'^ dtdrj^iuj) < I I \v{t,x,a)\^ dct,a{x)dt. (4.5) 
n(D) Jo Jo J D 

Then, because of (I3.15|) . it is easy to check that 77 := ?7„ (g) is a generalized incompressible 
flow, and moreover rj connects rj to 7. By integrating (|4.5p with respect to a, we obtain that 



5 Comparison of metrics and gap phenomena 

Throughout this section we shall assume that D = [0,1]*^. In [21], Shnirelman proved when 
d > 3 the following remarkable approximation theorem for Brenier's generalized (Lagrangian) 
flows: 

Theorem 5.1 If d > 3, then each generalized incompressible flow rj connecting i to h £ 
ST)iS{D) may be approximated together with the action by a sequence of smooth flows (^^(t, •)) 
connecting i to h. More precisely: 

(i) the measures r/i. := (5fc(-, x))#/i£) narrowly converge in il.{D) to ij; 

(ii) .s^xigk) = ■s/riVk) ^T{ri). 

This result yields, as a byproduct, the identity 

^(7i, Ih) = h) for all h G SDiff (D), d > 3. (5.1) 
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More generally the relaxed distance 5{r], 7) arising from the Lagrangian model can be compared, 
at least when rj = ji and the final condition 7 is induced by a map h G S{D), with the relaxation 
5* of the Arnold distance: 

:=inf jliminf5(i,/i„) : /i„ e SDiff(Z)), [ \hn - h\'^ dfio ^ o\ . (5.2) 

I n->oo Jjj J 

By ()3.7p and <\3.8h . we have > (^(7i,7h), and a gap phenomenon is said to occur if the 

inequality is strict. 

In the case d = 2, while examples of /i G SDiff (D) such that 6{i, h) = +00 are known |23j . 
the nature of (5* (h) and the possible occurrence of the gap phenomenon are not clear. 

In this section we prove the non-occurrence of the gap phenomenon when the final condition 
belongs to S{D), and even when it is a transport plan, still under the assumption d > 3. To 
this aim, we first extend the definition of 6^, by setting 

5*(7) := inf < liminf d{i, hn) ■ hn G SDiff(D), ~^ 7 narrowly > . (5.3) 

L n— >oo J 

This extends the previous definition (|5.2|) . taking into account that 7/^^ narrowly converge to 7/1 
if and only if /i„ ^ /i in L'^^fiu) (for instance, this is a simple consequence of [U Lemma 2.3]). 

Theorem 5.2 Ifd>3, then 5^(7) = ^(7^,7) for all 7 G r{D). 

The proof of the theorem, given at the end of this section, is a direct consequence of Theo- 
rem lS.ll and of the following approximation result of generalized incompressible flows by measure- 
preserving maps (possibly not smooth, or not injective), valid in any number of dimensions. 

Theorem 5.3 Let 7 G r(D). Then, for any probability measure rj on il.{D) such that 

{et)#r] = ^D VtG[0,r], (eo, er)#»7 = 7, 

and ^/riv) < there exists a sequence of flows {gk{t, ■))keN C W^'^ ([0, T]; L^(L')) such that: 

(i) gk{tr) G S{D) for all t G [0,r], hence rjf, := (^>gj#^z), with ^g^{x) = gk{-,x), are 
generalized incompressible flows; 

(ii) T]i^ narrowly converge in Q{D) to t] and s^T{gk) = -^TiVk) ~^ ^xiv)- 

Proof. The first three steps of the proof are more or less the same as in the proof of Shnirelman's 
approximation theorem (Theorem 15.11 in |24j ) . 

Step 1. Given e > small, consider the affine transformation of D into the concentric cube 

of size 1 — 4e: 

Te{x) := {2e,...,2e) + {l-Ae)x. 

This transformation induces a map from Q(D) into C([0, T]; D^) (which is indeed a bijection) 
given by 

f,{uj){t) ■.= T,{uj{t)) ycoeniD). 
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Then we define i)^ := {T^)jj,'q, and 

r], := (1 - AeYfi^ + r/o,„ 

where % £ is the "steady" flow in D\Di,: it consists of all the curves in D\D^ that do not move 
for < t < r. It is then not difficult to prove that rj^ ^ rj narrowly and a/xiVe) ~^ ■^Ti'n)^ 
e ^ 0. 

Therefore, by a diagonal argument, it suffices to prove our theorem for a measure t] which 
is steady near dD. More precisely we can assume that, if lo{0) is in the 2e-neighborhood of dD, 
then uj{t) = uj{0) for rj-a.e. uj. Moreover, arguing as in Step 1 of the proof of the above mentioned 
approximation theorem in [53], we can assume that the flow does not move for < t < e, that 
is, for rj-a.e. to, uj{t) = uj{0) for < t < e. 

Step 2. Let us now consider a family of independent random variables cji, 0^2, • • • defined in 
a common probability space (Z,Z,P), with values in C{[0,T], D) and having the same law rj. 
Recall that r] is steady near dD and for < t < e, so we can see uji as random variables with 
values in the subset of ^}(D) given by the curves which do not move for < t < e and in the 
2e-neighbourhood of the dD. By the law of large numbers, the random probability measures in 
n{D) 

1 ^ 

1=1 

narrowly converge to rj with probability 1. Moreover, always by the law of large numbers, also 

with probability 1. Thus, choosing properly we have approximated r} with measures i'tv 
concentrated on a finite number of trajectories uji{z){-) which are steady in [0,e] and close to 
dD. From now on (as typical in Probability theory) the parameter z will be tacitly understood. 

Step 3. Let (p G C^iW^) be a smooth radial convolution kernel with ¥?(x) = for > 1 
and ^p{x) > for |x| < 1. Given a finite number of trajectories uji, . . . ,lon as described is step 
2, we define 

a,(x) := (^ ^-^^M ^j dist(wi(0), dD) > e, 

«^(^) -^T.^ if dist(u;,(0), 51?) < e, 

where T is the discrete group of motions in generated by the reflections in the faces of D. 
It is easy to check that / = 1 and that supp(oj) is the intersection of D with the closed ball 
Be{oJi{0)). Define 

gi,t{x) := LJi{t) + {x - uJi{0)) Vz = 1, . . . , n. 

Let := (oi, • • • 1 o-N , 9i,t{x), . . . , gN,t{x)) and let us consider the generalized flow r}j~^ associ- 
ated to given by 

1 ^ 

/ /(^) ^^7V := T7 y] / ai{x)f{t^ gi,t{x))dx (5.4) 
Jn(D) ~1 Jd 



1=1 
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(that is, 77^ is the measure in the space of paths given by ai{x)5g. (^.^lix). The measure 

is well defined for the following reason: if dist(ti;j(0), 3D) < e we have gi^t{x) = x, and if 
dist{u;i{0), dD) > e and ai{x) > we still have that the curve t 1— > gi^t{x) is contained in D 
because ai{x) > implies jx — tijj(0)| < e and, by construction, dist{uJi{t),dD) > e for all times. 
Since the density p'^'^ induced by rjj^ is given by 

1 ^ 

i=l 

the flow ij^ is not measure preserving. However we are more or less in the same situation as in 
Step 3 in the proof of the approximation theorem in [24j (the only difference being that we do 
not impose any final data). Thus, by |24[ Lemma 1.2], with probability 1 

sup 1/0^ (t, re) — 1| ^ 0, 



sup p'\t, x)\ ^0 Va, 

x,t 

cT 



(5.5) 

dtp'\t,x)\' dtdx 



D Jo 



as — >■ 00. By the first two equations in (|5.5|) . we can left compose gi^t with a smooth correcting 
flow Ct^{x) as in Step 3 in the proof of the approximation theorem in |24j . in such a way that 
the flow fjj^ associated to := (ai, . . . , a^, C,^ o gi^[x), . . . o g^i(x)) via the formula 

analogous to (15. 4p is incompressible. Moreover, thanks to the third equation in (15. 5p and the 
convergence of £/t{i^n) to s^rivi)-! one can prove that ^^riVN) ~^ -^xiv) with probability 1. 

We observe that, since rj is steady for < t < e, the same holds by construction for t)^. 
Without loss of generality, we can therefore assume that Ct^ does not depend on t for t E [0,e]. 

Step 4. In order to conclude, we see that the only problem now is that the flow fj^v associated 
to is still non-deterministic, since if a; £ supp(aj) n supp(aj) for i 7^ j, then more that one 
curve starts from x. Let us partition D in the following way: 

D = DiLID2LI...LIDlLIE, 

where E is .if^-negligible, any set Dj is open, and all x G Dj belong to the interior of the 
supports of exactly M = M{j) < N sets aj, indexed by 1 < ii < • • • < < (therefore 
L < 2^). This decomposition is possible, as E is contained in the union of the boundaries of 
suppoj, which is ^'^-negligible. 

Fix one of the sets Dj and assume just for notational simplicity that ik = k for 1 < k < M . 
We are going to modify the flow f)^ in Dj , increasing a little bit its action (say, by an amount 
a > 0), in such a way that for each point in Dj only one curve starts from it. Given x € Dj, we 
know that M curves start from it, weighted with mass afc(x) > 0, and X]fe=i ^k{x) = 1. These 
curves coincide for < t < e (since nothing moves), and then separate. We want to partition 
Dj in M sets Ek, with 



^'^(Ek) = [ ttkix) dx, l<k<M 
Jd, 
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in such a way that, for any x € Ek, only one curve oo^ starts from it at time 0, G Dj for 

< t < £, and the map B x ^ ^xi^) ^ pushes forward ^'^\-Ej^ into a^^'^^Dj. Moreover, 
we want the incompressibihty condition to be preserved for all t G [0, e]. If this is possible, 
the proof will be concluded by gluing uj^ with the only curve starting from with weight 

The above construction can be achieved in the following way. First we write the interior 
of Dj, up to null measure sets, as a countable union of disjoints open cubes {Ci) with size 5i 
satisfying 

— E il^'^C.) < (5-6) 

i * 

with bi := min mina^. This is done just considering the union of the grids in given 

l<k<M Ci 

by for n G N, and taking initially our cubes in this family; if ()5.6I) does not hold, we 

keep splitting the cubes until it is satisfied (pi can only increase under this additional splitting, 
therefore a factor 4 is gained in each splitting). Once this partition is given, the idea is to move 
the mass within each Ci for < t < e. At least heuristically, one can imagine that in Ci the 
functions Ofc are almost constant and that the velocity of a generic path in Ci is at most of order 
5i/e. Thus, the total energy of the new incompressible fluid in the interval [0, e] will be of order 

V / r \i:ox{t)\''dtdx<-y^5fj^\c,) 

, Jcjo e , 

and the conclusion will follow by our choice of 6i. 

So, in order to make this argument rigorous, let us fix i and let us see how to construct our 
modified flow in the cube Ci for t £ [0,e]. Slicing Ci with respect to the first {d — l)-variables, 
we see that the transport problem can be solved in each slice. Specifically, if Ci is of the form 
x' + {0,5if, and we define 

:= / ak{x)dx, k = l,...,M, 
JCi 

whose sum is 5f, then the points which belong to Cf := + (0, 5iY~^ x have to move along 
curves in order to push forward ££'^i.C\ into a^^'^^Ci, where Jk are M consecutive open intervals 
in (0, 5j) with length 5\~'^m^. Moreover, this has to be done preserving the incompressibihty 
condition. 

If we write x = {x',Xd) G with x' = (xi, . . . ,Xd-i), we can transport the M uniform 
densities 

=^V(x* + {x'} X Jk) with x' G [0,Si]^'^-'^\ 

into the M densities 

afc(x',-)^V(x* + {x'} X [0,Si]) 

moving the curves only in the d-th. direction, i.e. keeping x' fixed. Thanks to Lemma 15.41 below 
and a scaling argument, we can do this construction paying at most M^b~'^6f/e in each slice of 
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Cj, and therefore with a total cost less than 

This concludes our construction. □ 

Lemma 5.4 Lei M > 1 be an integer and let 6i,...,6a/ : [0,1] — > (0,1] 6e continuous with 
= 1- Setting Ik = J^hkdt G (0,1], and denoting by Ji,...,Jm consecutive intervals 
of (0, 1) with length 1^, there exists a family of uniformly Lipschitz maps h{-,x), with h{t, •) G 
S{[0, 1]), such that 

h{l,-)#{XJ,^') = bk^\ k = l,...,M 

and 

^1 (/i) < , with b := min min6i. > 0. (5.7) 

^ ' - 62 ' l<k<M [0,1] ^ ^ 

Proof. We start with a preliminary remark: let J C (0, 1) be an interval with length / and 
assume that t ^ pj is a nonnegative Lipschitz map between [0,1] and L^(0, 1), with < 1 
and ptdx = I for all t G [0,1], and let f{t,-) be the unique (on J, up to countable sets) 
nondecreasing map pushing XJ-^^ to pt- Assume also that supppt is an interval and pt > r 
^^-a.e. on supp/Oj, with r > 0. Under this extra assumption, f{t,x) is uniquely determined for 
all x £ J, and implicitly characterized by the conditions 

/ pt{y)dy = ^\{0,x)n J), /(t, x) G supppj. 
Jo 

This implies, in particular, that f{-,x) is continuous for all x (z J. We are going to prove that 
this map is even Lipschitz continuous in [0, 1] and 

\^f{t,x)\ < ^^^EiEl for ^i-a.e. t G [0,1] (5.8) 
dt r 

for all X £ J. To prove this fact, we first notice that the endpoints of the interval supp pt (whose 

length is at least /) move at most with velocity Lip(p.)/r; then, we fix j; G J = [a, b] and consider 
separately the cases 

X £ dJ = {a, 6}, X G Int( J) = (a, 6). 
In the first case, since for any t G [0, 1] 

pt{y)dy = 0, / pt{y)dy = ^\j), 



JO 

and by assumption f{t,x) G supppf for any rr G J, we get supppf = [f{t,a),f{t,b)] for all 
t G [0, 1]. This, together with the fact that the endpoints of the interval supppj move at most 
with velocity Lip(p.)/r, implies ()5.8|) if x G dJ. In the second case we have 

pt{y)dye{0,.^\j)), 
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therefore f{t,x) £ Int(supppt) for all t € [0, 1]. It suffices now to find a Lipschitz estimate of 
\f{s,x) — f{t,x)\ when s, t are sufficiently close. Assume that f{s,x) < f{t,x): adding and 
subtracting Jq^^'^^ Pt{y) dy in the identity 



Pt{y)dy 



Ps{y)dy 



we obtain 



Pt{y)dy 



f{s,x) 



Ps{y) - pt{y)dy. 



Now, as f{s,x) belongs to supppt for \s — t\ sufficiently small, we get 

r\f{s,x)-f{t,x)\<Up{p.)\t 



s . 



connecting xJk-^^ 



This proves the Lipschitz continuity of f{-,x) and ([57 

Given this observation, to prove the lemma it suffices to find maps t 
to 6fc=Sf ^ satisfying: 

(i) supp/9j is an interval, and > miubk > b Jf^-a.e. on its support; 

[0)1] 

(ii) Up{p':) < ^ on [0, i], and Lip(p'=) < 2 on [i, 1]; 

M 

(ifi) J2 Pt for all t £ [0,1]. 
k=i 

Indeed, this would produce maps with time derivative bounded by (M — l)/(26) on [0, ^] and 
bounded by 2/6 on [i, 1], and this easily gives (j5.7p . 

The construction can be achieved in two steps. First, we connect xJk-^^ to Ik^^ in the 
time interval [0,^]; then, we connect Ik^^ to bk^^ in [^,1] by a linear interpolation. The 
Lipschitz constants of the second step are easily seen to be less than 2, so let us focus on the 
first interpolation. 

Let us first consider the case of two densities p^ = xJi and p^ = xjj, with Ji = (0, Zi) and 
J2 = (^1)0- 111 the time interval [0, r], we define the expanding intervals 



Jl,i = (0,/l + -/2), J'. 



2,t 



SO that Jfc,T = (0, 1) for A; = 1, 2, and then define 



pI-.- 



1 on{0,h-^h), 

h/l on (/i-iZi,/i + i/2) 

otherwise. 



pI 



By construction p^ > Ik on J^^t for k = 1, 2, pj + pi 

Lip(p'^) < 



{h--h,l), 

T 



(1 on{h + ^l2,l), 

_ otherwise. 
1, and it is easy to see that 



hh ^ I 
tI ~ At' 



(5.9) 
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We can now define the desired interpolation on [0, ^] for general M > 2. Let us define 

— — ! for i = 1,... ,M- 1, 

2(M-1) 

so that tM-i = ^- We will achieve our construction of on [0, ^] in M — 1 steps, where at each 
step we will progressively define on the time interval [ti^i,ti\. 

First, in the time interval [0, ti], we leave fixed pg •= XJ^-^^ for > 3 (if such k exist), while 
we apply the above construction in Ji U J2 to p^ and p"^. In this way, on [0, ii], := xJi-^^ is 
connected to ft\ := i^XJiUJ2-^^, and := XJ2-^^ is connected to p^^ := j^XJiUJa-^^- 

Now, as a second step, we want to connect p^^ to XJiUJ2UJ3-^^ for A; = 1, 2, 3, leaving 

the other densities fixed. To this aim, we define pj'^ := p^^ +p^j = XJiUJ2-^^- I^i the time interval 
[ti,i2]j we leave fixed Pq := xJk-^^ for /c > 4 (if such A; exist), and we apply again the above 
construction in Ji U J2 U J3 to pj^ and pj^ = XJ^-^^- In this way, on [ii,t2], Pti is connected 
to Pt2 ■= i^+lJli^ XJiUJzUJs-^^, and pf^ is connected to p| := i^+l+ijXJiUJaUJa-^^- Finally, it 
suffices to define pj := j^pf and p} := j^^pf ■ 

In the third step we leave fixed the densities pj^ for k > 5, and we do the same construction 
as before adding the first three densities (that is, in this case one defines p^^ ^ := pl^ + p^^ + Pt2 ~ 
XJiUJaUJa-^^)- In this way, we connect pj^^ to pj^^ := TT^^T^^fj^XJiUJaUJaUJ*-^^ and pf^ to 
Pts ■= X Jiu J2U J3U J4-^^ , and then we define p'l := j^rpf+ijpP^ for ^ = 1, 2, 3. 

Iterating this construction on [tj, tj+i] for i > 4, one obtains the desired maps t ^ p^. Indeed, 
by construction p\ > Ik on Jk^t, and Xlfcli = I- Moreover, by (15. 9p . it is simple to see that in 
each time interval one has the bound 

Lip(pn < . 

So the energy can be easily bounded by 1/6^ ( ^^^~^^" + l) ^ MV6^- □ 

Proof, (of Theorem 15. 2p By applying Theorem 15.31 to the optimal r) connecting i to 7, we can 
find maps gk G 5'(D) such that — > 7 narrowly and 

limsup%i,7gj < ^(7i,7)- 

Now, if d > 3 we can use ()3.12p . the triangle inequality, and the density of SDiff(L') in S{D) in 
the norm, to find maps G SDiff(L') such that 

Iimsup5(7i,7/,J < ^(7i,7) 
and 7/ij. — > 7 narrowly. This gives the thesis. □ 
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6 Necessary and sufficient optimality conditions 



In this section we study necessary and sufficient optimahty conditions for the generahzed geodesies; 
we shah work mainly with the Lagrangian model, but we will use the equivalent Eulerian- 
Lagrangian model to transfer regularity informations for the pressure field to the Lagrangian 
model. Without any loss of generality, we assume throughout this section that T = 1. 

The pressure field p can be identified, at least as a distribution (precisely, an element of the 
dual of ([0, 1] X D)), by the so-called dual least action principle introduced in [Hj. In order 
to describe it, let us build a natural class of first variations in the Lagrangian model: given a 
smooth vector field w{t, x), vanishing for t sufficiently close to and 1, we may define the maps 

: n{D) n{D) by 

S%uj,a){t) ■.= {e''"'uj{t),a), (6.1) 

where e^'^*x is the flow, in the (e, x) variables, generated by the autonomous field Wt{x) = w{t, x) 
(i.e. e*^^' = i and ^e^'^'x = w{t,e^'^^x)), and the perturbed generalized flows rj^ := {S^)^'q. 
Notice that 77^ is incompressible if divio^ = 0, and more generally the density p^^ satisfles for 
all times t G (0, 1) the continuity equation 

^p'^^it^x) + d\Y{wt{x)p'^^{t,x)) = 0. (6.2) 
ae 

This motivates the following definition. 

Definition 6.1 (Almost incompressible flows) We say that a probability measure u onVL{D) 
is a almost incompressible generalized flow if p^ G ([0, 1] x D) and 



-111 

5 - il|ci([0,l]xD) -• 



Now we provide a slightly simpler proof of the characterization given in |14j of the pressure 
field (the original proof therein involved a time discretization argument) . 

Theorem 6.2 For all rj, j £ T{D) there exists p G [C^([0, 1] x D)]* such that 

{P,p'' - l)(ci)*,ci < M(i^) -^^(r/,7) (6-3) 
for all almost incompressible flows f satisfying ()3.5p . 

Proof. Let us define the closed convex set C :={pe C\[0, 1] x D) : \\p - l\\ci < and the 
function : ([0, 1] x £>) ^ R+ U {+00} given by 

, ( inf {M(i/) : p" = p and ([33]) holds} if p G C; 

[ +00 otherwise. 

We observe that = 6 (?/,7). Moreover, it is a simple exercise to prove that (p is convex 
and lower semicontinuous in ([0, 1] x D). Let us now prove that has bounded (descending) 
slope at 1, i.e. 

hmsup-t-^-^ z^pU— < +00, 

p->l ||l-p||ci 
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By [m Proposition 2.1] we know that there exist < e < ^ and c > such that, for any p e C 
with II/? — 1 11(^1 < £, there is a Lipschitz family of diffeomorphisms gp{t,-) : D ^ D such that 

gp{t, ■)#PD = p{t, ■)liD, 

gp{t, ■) = Hot t = 0, 1, and the Lipschitz constant of {t, x) ^ gp{t, x) — x is bounded by c. Thus, 
adapting the construction in \14:\ Proposition 2.1] (made for probabihty measures in 0,{D), and 
not in Q{D)), for any incompressible flow r] connecting r] to 7, and any p £ C, we can define an 
almost incompressible flow v still connecting 77 to 7 such that p"^ = p, and 

^iH < ^liv) + c'Wp - l||ci(l + ■f^iiv)), 

where d depends only on c (for instance, we define v := G^r}, where G : Q{D) Q{D) is the 

map induced by gp via the formula (w(t),a) 1-^ {gp(t,u;{t)), a)). In particular, considering an 
optimal ?7, we get 

0(p) < ^{1) + c\\p - l\\ci (1 + f{v, 7)) (6.4) 

for any p £ C with ||p — l||c7i < £■ This fact implies that <j) is bounded on a neighbourhood of 1 
in G. Now, it is a standard fact of convex analysis that a convex function bounded on a convex 
set is locally Lipschitz on that set. This provides the bounded slope property. By a simple 
application of the Hahn-Banach theorem (see for instance Proposition 1.4.4 in [3]), it follows 
that the sub differential of at 1 is not empty, that is, there exists p in the dual of G^ such that 

{p, p - i)(ci)*,ci < Hp) - <?^(i)- 

This is indeed equivalent to (j6.3p . □ 
This result tells us that, if 77 is an optimal incompressible generalized flow connecting 77 to 7 

—2 

(i.e. ^1(77) = 6 (77,7)), and if we consider the augmented action 

ii^)-= i [ h^{t)\^dtdu{uj,a)-{p,p^-l), (6.5) 
Jn{D) Jo ^ 

then r\ minimizes the new action among all almost incompressible flows v between r/ and 7. 
Then, using the identities 



dedt ^ ' 



= ^w{t, uj{t)) = dtw{t, uj{t)) + V^w{t, uj{t)) ■ ijj{t) 

e=0 



and the convergence in the sense of distributions (ensured by (|6.2p ) of (p^^ — l)/e to —divw as 
e J, 0, we obtain 



= / / uj{t) ■ ■^w{t,uj{t)) dtdr]{u!,a) + {p,divw). (6.6) 
=0 Jn{D) Jo dt 



As noticed in [TJj , this equation identifies uniquely the pressure field p (as a distribution) up to 
trivial modifications, i.e. additive perturbations depending on time only. 
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In the Eulerian-Lagrangian model, instead, the pressure field is defined (see (2.20) in |16j ) 
and uniquely determined, still up to trivial modifications, by 

Vp{t,x) = -dt (^j v{t,x,a)dct,x{a)^ - div v{t,x,a) ^ v{t,x,a) dct,x{a)^ , (6.7) 

all derivatives being understood in the sense of distributions in (0, 1) x D (here (c, v) is any 
optimal pair for the Eulerian-Lagrangian model). We used the same letter p to denote the 
pressure field in the two models: indeed, we have seen in the proof of Theorem 14.11 that . writing 
V = Va f^D, the correspondence 

V ^ ict,a^v2a) with cj^„ := iet)#ri^, v^^^c^^^ := {et)#[^[t)j]^) 

maps optimal solutions for the first problem into optimal solutions for the second one. Since 
under this correspondence ()6.7p reduces to ()6.6p . the two pressure fields coincide. 

The following crucial regularity result for the pressure field has been obtained in [S], im- 
proving in the time variable the regularity dxj) G A^ioc((0, 1) x D) obtained by Brenier in 

m- 

Theorem 6.3 (Regularity of pressure) Let (c, v') he an optimal pair for the Eulerian-La- 
grangian model, and let p he the pressure field identified by (HZD- Thendx.p G -L^ ((0' 
and 

peLl^{{Q,l)-BVUD)) C LL((0,l);Lf/^('"')(I))). 

Ln the case D = T'^ the same properties hold glohally in space, i.e. replacing BV\oc{D) with 
BV{T'^) and L'11^'^~^\d) with L'^/^'^-^\T'^). 

The Ljq^ integrability of p allows much stronger variations in the Lagrangian model, that 
give rise to possibly nonsmooth densities, which may even vanish. 

From now one we shall confine our discussion to the case of the fiat torus T'^, as our arguments 
involve some global smoothing that becomes more technical, and needs to be carefully checked 
in more general situations. We also set /if = ^^jd and denote by dj the Riemannian distance in 
T"^ (i.e. the distance modulo 1 in M'^/Z'^). In the next theorem we consider generalized fiows u 
with hounded compression, defined by the property p'^ G L°° ((0, 1) x D). 

Theorem 6.4 Let r] be an optimal incompressible flow in T'^ between rj and 7. Then 

{p, - 1) < s^iii^) - £^iiv) (6.8) 
for any generalized flow with bounded compression v between r] and 7 such that 

p^{t, •) = 1 for t sufficiently close to 0, 1. (6-9) 
If p & -^'^([0, 1] X T'^), the condition ()6.9p is not required for the validity of (16. Sp . 
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Proof. Let J := •) 7^ 1} (0, 1) and let us first assume that p"^ is smooth. If — 111(^1 < 

1/2, then the result follows by Theorem 16.21 If not, for e > small enough (1 — e)ri + ei^ is a 
slightly compressible generalized flow in the sense of Definition 16.11 Thus, we have 

e{p, - 1) = {p, - 1) < i^i((l - s)^ + ei^) - M(^) = e {s^i{u) - ^(^7)) , 

and this proves the statement whenever is smooth. 

If p^ is not smooth, we need a regularization argument. Let us assume first that p^ is smooth 
in time, uniformly with respect to x, but not in space. We fix a cut-off function x £ ^^{Q, 1) 
strictly positive on a neighbourhood of J and define, for y G W^, the maps T^^y : (}{T'^) 0(T'^) 
by 

Te,y{uj, a) := {to + eyx, a), {cu, a) G J^(T'^). 

Then, we set := f^a{T^,y)#i^(t>{y) dy, where (j) : Mf^ —)■ [0, +00) is a standard convolution 
kernel. It is easy to check that still connects rj to 7, and that 

p'^^(t,-)=/o"(t,-)*4x(i) VtG [0,1], 
where (pe{x) = e~'^(t){x / e). Since 

lim ^2/1(1/;.) = lim / / / \Lo{t) + £yx{t)\'^ dtdu{u),a)(j){y) dy = S2ti{v) 
e-lO JK.'i Jn{T'i) Jo 

we can pass to the limit in ()6.8p with in place of v, which are smooth. 

In the general case we fix a convolution kernel with compact support (p{t) and, with the 
same choice of x done before, we define the maps 

T,{co, a)it) := ( / Lo{t - sex{t)Ms) ds, a). 
Jo 

Setting i/'^ = (T^)^i', it is easy to check that ■s^xiy') and that 

p''-{t,x)= [ p''{t-sex{t),x)ip{s)ds 
Jo 

are smooth in time, uniformly in x. So, by applying ()6.8p with in place of v, we obtain the 
inequality in the limit. 

Finally, if p is globally integrable, we can approximate any generalized flow with bounded 
compression v between r] and 7 by transforming co into cooips^ where : [0, 1] [0, 1] is deflned 
by ipe{t) := j^2^ Jo X[e,i-e]{^)ds (so that -i/'e is constant for t close to and 1). Passing to the 
limit as e I we obtain the inequality even without the condition p'^{t, •) = 1 for t close to 0, 1. 

□ 

Remark 6.5 (Smoothing of flows and plans) Notice that the same smoothing argument 
can be used to prove this statement: given a flow 77 between i] = rja pr and 7 = 7a /^t 
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(not necessarily with bounded compression), we can find flows with bounded compression rj^ 
connecting r]^ := {rja) * (j)s fJ-j to := (7^) * c^e ® fij, with s^xiv^) = -^Tiv) 

[ [ rs{T,uj)dTdri{uj,a) = f [ r{T,u)) dr drj" {uj,a) Vr G ([0, 1] x T'') 

(where, as usual, r^{t,x) = r{t, •) * (j)fr{x)). In order to have these properties, it suffices to define 

where C7^(a;,a) = (w + z,a). Notice also that the "mollified plans" r]'^ , 7^ converge to 77, 7 
in (r(T'^),5): if we consider the map 5^ : T'^ — > Q{T'^) given by x 1— > ujx(t) := x + ety, the 
generalized incompressible flow i/^ = ® /iT, with 

{S'y)#\a4>{y)dy, 

connects in [0, 1] the plan A = Aa (8> /^t to = (Aa * (^e) (8> /^t, with an action equal to 

e^Rd \y?<i>ky)dy- 

In order to state necessary and sufficient optimality conditions at the level of single fluid 
paths, we have to take into account that the pressure field is not pointwise defined, and to 
choose a particular representative in its equivalence class, modulo negligible sets in spacetime. 
Henceforth, we define 

j5(i, x) := liminf ^^(t, x), (6.10) 
ej,0 

where, thinking of p(t, •) as a 1-periodic function in M*^, is defined by 

Veit, x) := (27r)-'^/2 / p{t, x + ey)e~\y\"/^ dy. 



Notice that Ps is smooth and still 1-periodic. The choice of the heat kernel here is convenient, 
because of the semigroup property p^+g' = (Pe)e'. Recall that p is a representative, because at 
any Lebesgue point x of p(t, •) the limit oi p^{t,x) exists, and coincides with p{t,x). 

In order to handle passages to limits, we need also uniform pointwise bounds on p^; therefore 
we define 

M/(x) := sup (27r)-'^/2 / \f\{x + ey)e-\y\'/^dy, f e L^T''). (6.11) 

e>0 JRd 

We will use the following facts: first, 

Mf, = sup \fe\e' < SUp(|/|,),. < SUp |/|, = Mf 
e'>0 e'>0 r>0 

because of the semigroup property; second, standard maximal inequalities imply HM/H^p^jd) < 
CpII/||lp(T'*) fo'^ s-ll P > 1- Setting Mp{t,x) := Mp{t, •){x), by Theorem 16.31 we infer that Mp G 
Lf^^((0, l),L''/'^-^(T"')), so that in particular Mp G ^[^^^((0, 1) x T'^). This is the integrability 
assumption on p that will play a role in the rest of this section. 
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Definition 6.6 (g-minimizing path) Letu£H'^{{0,l)-D) with Mq{T,uj) £ L^{0,1) 
We say that lo is a g-minimizing path if 

-\^{t)\^ -q{T,uj)dT < / -\^{t)+5{t)\^ -q{T,UJ + 5)dT 



for all 5 € {{0,1); D) with Mq{T,uJ + 6) G L^{0,1). 
Analogously, we say that oj is a locally (/-minimizing path if 

£ i|d;(r)|2 - q{T, uj) dr < £ ]^\Co{t) + 5{t)\'' - q{T, u: + 5) dr (6.12) 

for all [s,t] C (0, 1) and all 5 £ {{s,t);D) with Mq{T,uj + 6) £ L^{s,t). 

Remark 6.7 We notice that, for incompressible flows rj, the (resp. L^^^) integrability of 
Mq{T, Lo) imposed on the curves uj (and on their pertm'bations u + 6) is satisfied ry-a.e. if 
Mq G ((0, 1) X T'^)) (resp. Mq G LJ^^ ((0, 1) x T"^)); this can simply be obtained first noticing 
that the incompressibility of t] and Fubini's theorem give 

/ / f{T,uj)dTdri{u;,a) = / / f{T,x)dfijd{x)dT 

for all nonnegative Borel functions / and all intervals J C (0, 1), and then applying this identity 
to / = Mq. 

Theorem 6.8 (First necessary condition) Let ij — Va ^ f^'^ ^''^V optimal incompressible 
flow on T"^. Then, rj is concentrated on locally p-minimizing paths, where p is the precise 
representative of the pressure field p, and on p-minimizing paths if Mp G -Z>^([0, 1] x T*^). 

Proof. With no loss of generality we identify T"^ with M'^/Z'^. Let r] be an optimal incompressible 
flow and [s,t] C (0, 1). We fix a nonnegative function x G ^^{0, 1) with {x > 0} = {s,t). Given 
6 G {[s,t];T'^), y G M"^ and a Borel set E C n{T'^), we define T,^y : 0(T'^) ^ ^{T"^) by 



Te,y{uJ,a) :-- 



{lo, a) if w ^ E; 

{lo + 5 + eyx, a) if lo £ E 



(of course, the sum is understood modulo 1) and Ve^y := {T^^y)^!^. 

It is easy to see that Vf^^y is a flow with bounded compression, since for all times r the curves 
uj{t) are either left unchanged, or translated by the constant 6{t) +£yx{T)^ so that the density 
produced by Ugr^y is at most 2, and equal to 1 outside the interval [s,t]. 

Therefore, by Theorem 16.41 we get 

/ / p{p'^^-'' - I) dr d^j < I £/i{lo + 6 + eyx) - J!^i{Lo)dT]{LO,a). 
Rearranging terms, we get 



-p{T,ijj)dTdr]{uj,a) < ]^\uj + 5 + eyx\^ - p{t , uo + 5 + eyx) dr 



dr]{uj, a). 
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We can now average the above inequality using the heat kernel cf>{y) = (27r) "^/^e 1^1^/^, and we 
obtain 

t I 

-|u;p — p{t,uj) dr dr]{uj,a) 

~ Ie ^I^"^^"^^^^'^^'^'^*^^''^^" / ^^'cM*^'^''^"^'^''^'^ dr]{uj,a). 

Now, let V C Hq ([s,t];T'^) be a countable dense subset; by the arbitrariness of E and Re- 
mark [^3] we infer the existence of a rj-negligible Borel set B C Cl(T'^) such that Mp{T,uj) € 
L^{s, t) and 

^\uj\^ -p{T,uj)dT < ff ^\uj + 6 + eyxl"^ dT(f>{y) dy - I p^^(^^){t,uj + 6) dr 

J K.*^ J s J s 

holds for all e = 1/n, 5 G P and {u},a) G fi(T'^) \ By a density argument, we see that the 
same inequality holds for all e = 1/n, 6 G Hi ([s,t];T'^), and oj G fi(T'^) \ 

Now, \iMp{T,uj + 5) G since 6 £ Hi [[s,t\;T'^) we have that Mp{T,uj + 5) G 

and we can use the bound \pfr\ < Mp to pass to the limit as e | to obtain that (|6.12|) holds 
with q = p. 

The proof of the global minimality property in the case when p G L^([0,T] x T*^) is similar, 
just letting 6 vary in Hi ([0, 1]; T'^) and using a fixed function x G ^^([0, 1]) with x(0) = = 
and X > in (0, 1). □ 

In order to state the second necessary optimality condition fulfilled by minimizers, we need 
some preliminary definition. Let q £ {[^^A ^ ^) us define the cost Cg* : D x D ^ 

of the minimal connection in [s, t] between x and y, namely 

:=inf 1^ ^\u;{t)\^ - q{T,uj) dr : uj{s) = x, uj{t) = y, Mq{T,uj) e 1^8,1)^ (6.13) 

with the convention Cq^(x,y) = +oo if no admissible curve uj exists. Using this cost function 
Cg'*, we can consider the induced optimal transport problem, namely 

W^,M(/xi,/i2) :=inf I / cl'\x,y)dX{x,y) : A G r(^i, /xa), (cf )+ G L1(A)| , (6.14) 
' IJdxd J 

where r(/ii,/i2) is the family of all probability measures A in D x D whose first and second 
marginals are respectively /ii and fi2- Again, we set by convention W_s,t(/xi, //2) = +oo if no 

Cq 

admissible A exists. 

Unlike most classical situations (see |26j). existence of an optimal A is not guaranteed because 
Cg* are not lower semicontinuous in D x L>, and also it seems difficult to get lower bounds on 
Co*. It will be useful, however, the following upper bound on W s,t: 

Cq 

Lemma 6.9 If Mq G L^([s,t] x T'^) there exists a nonnegative fif-integrable function Kq^ 
satisfying 

cf{x, y) < Kl^\x) + ^^\y) Vx, y G T^. (6.15) 
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Remark 6.10 By (l6J5]l we deduce that, if Kg'* G L^{fii + ^2), then {c/)+ G L^{X) for all 
A G r(/ii, /X2) and we have 



c'/{x,y)dX{x,y) < K'/{w) d{fii + fi2){w) VA G r(/ii, //s). 
In particular, s,t(/ii, /U2) as defined in (I6.14p is not equal to +00. 

Proof. Assume s = and let I = t/2. Let us fix x, y G T*^; given 2; G T'^ we consider the 
projection on T"' of the Euclidean path 



x + j{z-x) ifTG[0,/]; 
z+^{y-z) ifTG[/,t]. 



This path leads to the estimate 



y) < dU^^')+dUz,y) ^ £ ^^^^^^ X + - x)) dr + J' Mq{r, z+^-^{y-z)) dr. 

Tpci ^ Td 



By integrating the free variable z with respect to //t, since < on T x T , we get 



c°'*(x,y)< ^ + 
Therefore, the function 



f [ Mq{T,x + j{z-x)) + Mq{l + T,z + j{y-z))dTdfiT:{z) 
JT'^ Jo ' ' 



-^a'*H:= ^+/" [ Mq{T,w + ^{z-w))+Mq{l + T,z+^{w-z))dTdfiT{z) (6.16) 
4/ Jjd Jq I I 

fulfils ()6.15p . It is easy to check, using Fubini's theorem, that Kq'^ is ^T-integrable in T*^. Indeed, 
/ Kl^\w)diJLr{w) = ^+ f [ [ Mq{T,w + j{z-w))dTdnj{z)dfiT{w) 
+ 111 Mq{l + T,z + j{w - z))dTdfij{w)dfij{z) 

Jjd Jjd Jq L 

= JJ+[ [ I Mq{T,w + jy)dTdnT{y)dfij{w) 

4t Jjd Jjd Jq I 

+ Mq{l + T,z + jy)dTdfij{z)dfij{y) 

Jjd Jjd Jq I 



d 

4/ + 



mMq{T, w + jy) + Mq{l + t,w + jy) dnj{w) diJ,j{y) dr 



+ / / Mq{T,w) dfij{w) dr < +00. 
'0 -'T'' 



□ 
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In the proof of the next theorem we are going to use the measurable selection theorem (see fL9\ 
Theorems III. 22 and III. 23]): if {A, A, v) is a measure space, X is a Pohsh space and E d Ax X 
is Au 5(X)-measurable, where Au is the i/-completion of A, then: 

(i) the projection tta{E) of E on A is „4^-measurable; 

(ii) there exists a {Au, ;B(X))-measurable map a : vr(£') — > X such that (x, cr(x)) G for ly-a.e. 
X e 7rA(^). 

The next theorem wiU provide a new necessary optimahty condition involving not only 
the path that should be fohowed between x and y (which, as we proved, should minimize the 
Lagrangian Cp in p.8|) ). but also the "weights" given to the paths. We observe that, when 
a variation of these weights is performed, new flows r) between r/ and 7 are built which need 
not be of bounded compression, for which {et)^f} might be even singular with respect to fij; 
therefore we can't use directly them in the variational principle ()6.8p : however, this difficulty 
can be overcome by the smoothing procedure in Remark [ 



Theorem 6.11 (Second necessary condition) Letr] = Ty^^/ix be an optimal incompressible 
flow on T'^ between rj and 7. Then, for all intervals [s,t] C (0, 1), W s_,t [r]a,"/a) £ IR and the plan 



{es,et)#r]^ is optimal, relative to the cost Cp defined in ()6.13p . for fij-a.e. a. 
Proof. Let [s,t] C (0,1) be fixed. Since 

/ Cp^\x,y)d{es,et)#r]^dfij{a) < [ [ [ l-\dj{T)\^ - p{t,uj) drdrj^iuj) dfij{a) 



= it-s)6\7j,j), (6.17) 

it suffices to show that 

{t-s)fiv,7)< [ H^,M(r/^,7a)rf^T(a). (6.18) 

We are going to prove this fact by a smoothing argument. We set r/* = r]^(d) fij, 7* = 7* ® /ii, 
with 77^ = {es)#r]^, 7* = {et)#r]^. Recall that Remark l3 . 1 1 gives 

6iri,ri') = s5(?7,7), 5(7*,7) = (1 -*)%,7)- 
First, we notice that Lemma 16.91 gives 



Jd -\p\ Jfd Jjd 



(6.19) 



jd 



K^^-^{uj)dfij{uj) < +00. 



We also remark that, since r 1— > ||pe(T, •)||oo is integrable in (s,t), for any e > the cost Cp'^ is 
bounded both from above and below. Next, we show that 

Cp'\x,y) > limsupc^f(x,?/) V(x,y) G T"^ x T"'. (6.20) 

ej,0 
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Indeed, let oj G ([s,t];T^) with u}{s) = x, uj{t) = y and Mp(r, w) G L^{s,t) (if there is no 
such uj, there is nothing to prove). By the pointwise bound < Mp and Lebesgue's theorem, 
we get 

1 /■* 1 

-\uj{t)\'^ - p{T,u;)dT = lira / -|u;(r)|2 - ^^(t, w) dr. 

By the L^(L°°) bound on Mp^, the curve ijj is admissible also for the variational problem defining 
Cp^*, therefore the above limit provides an upper bound on limsup^ Cp^{x, y). By minimizing with 
respect to to we obtain (|6.2U|) . 

By (16. 191) and the pointwise bound p > —\p\ we infer that the positive part of W 7a) is 

/ix-integrable. Let now 5 > be fixed, and let us consider the compact space X := i^(T'^ x T*^) 
and the ^(T*^)^^ ® ;B(X)-measurable set 

i^:= |(a. A) GT'^xX: AGr(r?^,7*), / cf (x, y) dA < 5 + (l^^s^r?^, 7*) V 1 

(we skip the proof of the measurability, that is based on tedious but routine arguments) . Since 
W a,t{r]'i,-^V) < +00 for fij-a.e. a, we obtain that for fij-a.e. a G T'^ there exists A G r(r/* 7*) 

with (a, A) G E. Thanks to the measurable selection theorem we can select a Borel family 
a Aa G ^{T'^ X T'^) such that Aa G r(??^,7*) and 



for fij-a.e. a G T'^. 



By Lemma 16.91 and Remark 16.101 we get 

(x, y) < K;'^\x) + K^'^iy) < K;\x) + K;\y) Vx, y G T'^ 

and 

/ / K;^\x) + K;4y) dXa dfij{a) =11 K;^' d{r,i + 7*) dfij{a) < +00 

Jjd Jjdy.jd Jjd Jjd 

(we used the pointwise bound Mp/r < Mp and the fact that q Kq'^ has a monotone dependence 
upon Mq, see ()6.16p ). Therefore (I6.20p and Fatou's lemma give 

6+1 (vF^s,t(7?^,7*) V-^) (i/iT(a) > limsup / / ct^!{x,y)dXadfij{a). (6.21) 

Still thanks to the measurable selection theorem, we can find a Borel map {x,y,a) 1— > uja'J G 
C {[s, t];T'i) with <|(s) = X, u^alUt) = y, Mp,{T, ujI^^) G L\s, t) and 

/■* 1 

/ 2^^a,e\'^ - Pe{T,^a,e) dr < 5 + c^f (x, y) for Aa nj-a.e. {x, y, a). 

J S 
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Let = (8) be the push-forward, under the map (x, y, a) i— > Wa'l , of the measure Xa ® //t; 
by construction this measure fulfils (es)^A^ = rj^, (et)#A^ = 7*, (because the marginals of Xa 
are rj^ and 7*), therefore it connects 77^ to 7* in [s,t]. Then, from ()6.2ip we get 

2S+ [ (w^s.t{rji,ji)V -^)dnjia)>\imsvLp [ [ J|c0(r)|2 - p,(t, dr dA"(cu, a). 

Eventually, Remark 16.51 provides us with a flow with bounded compression A^ connecting rf'^ 
to 7*'^ in [s, t] with 

2(5+ / ('w .,t(?7a,7a) V (i;UT(a) > limsup / / ^|u;(r)|^ - p(t, cj) dr ^^"(0;, a). 

(6.22) 

Since rf''^ — > 77* and 7*'^ — > 7* in (r(T'^),5), we can find (by scaling r\ from [0,s] to [0,Se] and 
from [i, 1] to [tg, 1], and using repeatedly the concatenation, see Remark |3.2|) generalized flows 
v"^ between 7 and i] in [0, 1], t s, | t satisfying: 

(a) connects r\ to jf in [0, s^], r/** to rf'^ in [sg, s], 7*'^ to 7* in [t, tg], 7* to 7 in [t^, 1] and is 
incompressible in all these time intervals; 

(b) the restriction of to [s,t] coincides with A"^; 

(c) the action of v'' in [0,,s] converges to b [i],!]^) = (77,7), and the action of i/^ in [t, 1] 
converges to 5 (7^,7) = (1 — t)'^6 {'i],'y)- 

Since is a flow with bounded compression connecting r/ to 7 we use (|6.8p and the incompress- 
ibility in [0, 1] \ [s,t] to obtain 

/ / l\djiT)\^ dTdu%u;,a) - [ [ P{t,u;) dr du'{uj,a) >f{r],-f) (6.23) 

for all e > 0. Taking into account that (b) and (c) imply 

Q(T<*) Jo ^ 

and ^ 

/ / l|c:;(T)|2dTdi.^(a;,a)^(l-t)5'(r/,7), 

from and we get 

25 + ^^(W^A<.1\) V -^) d^T(a) > (1 - s - (1 - t))^'(7?,7) = (t - 5)5' (r?, 7). (6.24) 
Letting 5 | we obtain the //T-integrability of s,t (r;^, 7*) and (|6.18p . □ 
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A byproduct of the above proof is that equahties hold in (I6.17P , (I6.18P , and therefore 

T'i Jn{T'i) \Js 2 



(^{t)]^ - P{r,uj) dr - c/{uj{s),uj{t)) ] dr]^{uj) dnj{a) (6.25) 



= / / [ h'^i^)]^ -Pi'r,'^)dTdr]^{uj)dfxj{a)- [ W^sXif^,-/^) dnj{a) = 0. 

This yields in particular also the first optimality condition. However, as the proof of Theo- 
rem [6TTT] is much more technical than the one presented in Theorem 16.81 we decided to present 
both. 

Now we show that the optimality conditions in Theorems 16.81 and 16.111 are also sufficient, 
even in the case of a general compact manifold without boundary D. 

Theorem 6.12 (Sufficient condition) Assume that rj = rj^® is a generalized incompress- 
ible flow in D between rj and 7, and assume that for some map q the following properties hold: 

(a) Mq £ L^((0, 1) x D) and Tj is concentrated on q-minimizing paths; 

(b) the plan {eo,ei)^r]^ is optimal, relative to the cost Cg'^ defined in (I6.13p . for ^D-a.e. a. 
Then rj is optimal and q is the pressure field. In addition, if (a), (b) are replaced by 

(a') Mq £ Ll^iiO, 1) x D) and rj is concentrated on locally q-minimizing paths; 

(b ') for all intervals [s, t] C (0, 1), the plan (e^, et)^r]g^ is optimal, relative to the cost Cg* defined 
in ()6.13p . for fiD-a.e. a, 

the same conclusions hold. 

Proof. Assume first that (a) and (b) hold, and assume without loss of generality that 
q{t, ■) dfi£) = for almost all t S (0,1). Recalling that, thanks to the global integrability 
of Mq, any generalized incompressible flow v = between r/ and 7 is concentrated on 

curves uj with Mq{T,uj) G i^^(0, 1) (see Remark 16. 7p . we have 

J^i{u) =11 I hdj]'^ - q{T,oj)dTdVa{(^)dfiD{a) (6.26) 
Jd Jn{D) -Jo ^ 

> c^g''^{x,y)d{eo,ei)#UadfiD{a) > / 0,1 (r?a, 7a) c^/^d(«)- 

When v = rj the flrst inequality is an equality, because rj is concentrated on g-minimizing paths, 
as well as the second inequality, because of the optimality of the plan {eQ,ei)^ri^. This proves 
that 77 is optimal. Moreover, by using the inequality in ()6.26p with a flow v with bounded 
compression, one obtains 

£/i{u)>£/i{rj) + {q,p''-l). 

Considering almost incompressible flows u arising by a smooth perturbation of rj as described 
at the beginning of this section (see (16. ip in particular), the same argument used to obtain (16. 6p 
gives that q satisfies (|6.6p . so that q is the pressure field. 
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In the case when (a)' and (b)' hold, by locahzing in all intervals [s,t] C (0,1) the previous 
argument (see Remark I3.2p . one obtains that 



(t-s) [ [ ^|wp(iT(iT7(w,a) = (5^(7^,7t), 
Jn(D) Js ^ 



in{D) Js 

where 7^ = (e^, vr£))^?7 and 74 = {et,'n'D)#'n- Letting s [0 and t j 1 we obtain the optimality of 
rj. □ 

A byproduct of the previous result is a new variational principle satisfied, at least locally in 
time, by the pressure field. Up to a restriction to a smaller time interval we shall assume that 
Mp E ([0,1] X T"'). 

Corollary 6.13 (Variational characterization of the pressure) Let 7 G r(T'^) and let 
p be the unique pressure field induced by the constant speed geodesies in [0, 1] between r] = 
Va fJ'J and 7 = 7a ® /^T- Assume that Mp € ([0, 1] x T*^) and, with no loss of generality, 
jjd p{t, •) dfij = 0. Then p maximizes the functional 



Jjd 1 Jq Jjd 



(r, x) dfij{x) dr 

among all functions q : [0,1] x T'^ ^ R with Mq e L^{[0, 1] x T"*). 

Proof. We first remark that the functional ^ is invariant under sum of functions depending on 
t only, so we can assume that the spatial means of any function q vanish. 
From ()6.25p we obtain that 



/ 

Jt 



W^o,i{r]a,ja)dfij{a) = I I I ^\uj{T)f -p{T,uj)dTdr]^{uj)dnj{a). 



T'* Jf^ in(T'*) 



By the incompressibility constraint, in the right hand side p can be replaced by any function 
q whose spatial means vanish and, if Mq £ L} ([0, 1] x T'^), the resulting integral bounds from 
above /^^ W^o,i (r/a, 7a) diif{a), as we proved in ()6.26p . □ 
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